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Z boson transverse momentum 
spectrum from the lepton angu- 
lar distributions 
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Abstract — In view of recent discussions concerning the possibly limiting energy resolution systematics 
on the measurement of the Z boson transverse momentum distribution at hadron colliders, we propose a 
novel measurement method based on the angular distributions of the decay leptons. We also introduce a 
phenomenological parametrization of the transverse momentum distribution that adapts well to all currently 
available predictions, a useful tool to quantify their differences. 
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1 Introduction 



The transverse momentum distribution of heavy particles at hadron colliders is a longstanding subject, first 
discussed in the context of QCD. Formalisms were developed to predict this distribution, based on analytical 
or numerical methods (soft gluon resummation theory |Q]- 4fl, or parton shower Monte Carlo programs J3,@] 
respectively). In both cases, the shape of the distribution is predicted qualitatively, but the full result depends 
on a limited number of free parameters which need to be extracted from measurement. 



The distribution is of physical interest for many reasons. Firstly, the measurement of the overall vector bo- 
son production cross sections is considered an important test of perturbative QCD, theoretical predictions 
now being available up to NNLO |7J]. In the context of the total cross section measurement, the kinematic 
cuts imposed on the decay leptons, reflecting the detector geometric acceptance, require that the observed 
event rate be corrected by a factor compensating this loss of acceptance. The fraction of lost events must be 
precisely controled so that the final result contains no significant bias. This in turn implies that the lepton 
kinematic distributions, and hence the vector boson ones from which they derive, need to be known both 
inside and outside the selected region. 



Another application is the precise measurement of the W boson mass QSQ . The decay lepton transverse mo- 
mentum distribution, or the W boson transverse mass distribution from which this fundamental parameter 
is extracted, is a complicated quantity resulting in part from the W boson transverse momentum spectrum, 
dow/dpr- The increasing precision of the measurements of My? puts ever stronger constraints on the knowl- 
edge of dow/dpr- An important tool for constraining this distribution is the study of the Z boson transverse 
momentum spectrum, doz/dpr- 



The recent measurements performed at the Tevatron are increasingly sensitive to the detector energy res- 
olution, which needs to be precisely "subtracted", or unfolded from the observed distribution to derive an 
estimate of the true one. This issue will become much more important at the LHC, given the high expected 
statistics. An alternative variable, a T , was introduced recently [9] as a replacement for p T . Its advantage 
is that it is negligibly sensitive to the energy resolution, while still a good probe of resummation or parton 
shower mechanisms. 



In this paper, we propose a novel method that shares the insensitivity of aj to the energy resolution, while 
remaining a true pj measurement. The method is based on the measured angular distributions of the Z boson 
decay leptons, which, together with the well known Z boson mass, are sufficient to extract the pj distribution. 



In the following, we first introduce a convenient parametrization of the p T distribution. It depends on three 
intuitive parameters, and adapts well to the available predictions. It represents a practical tool to quantify 
differences between predictions, as well as for the measurement itself. We then outline the measurement 
method, and give examples of its performance in a simplified form. We conclude with some caveats and 
perspectives concerning the use of this method in future measurements. 
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2 A parametrized form for the p T distribution 

The parametrization we propose relies on a number of simple arguments. Consider Z boson production at 
high energy, and at given mass and rapidity, so that the partem momentum fractions at the hard vertex are 
small and fixed. In the low transverse momentum region, the repeated gluon emission in the initial state 
generates a gaussian transverse momentum distribution. Along both the x and y axes, this "random walk" 
leads to a distribution proportional to 
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f{p x y,o PS ) dp XJ ~ e 2a es dp. 



The Ops parameter represents the spread of the px.y distribution after all emissions and, in a naive picture, 
could be seen as representing the average number of emitted gluons times their average transverse momen- 
tum: Ops ~ \/Ng x Pxy- Moving to polar coordinates, the distribution becomes: 



fiPxl o>s) f{j>y, Ops) dp x dp y ~ e 2<T « e 2a « dp x dp y 

p T dp T d<$> = g\(p T ;Ops) dp T - 



2 
Pf 



after a trivial azimuthal integral. At higher pj, the shape is dominated by a power law behaviour representing 
the parton density functions (PDFs) and the perturbative matrix element: 

8l(PT;a) ~ l/p£. 

The transition between the two descriptions is controled by a parameter n, defined such that pj atch = nx Ops. 
As for the definition of the Crystal Ball function ifioh . it satisfies smoothness conditions (the function and its 
derivative are continuous). The complete parametrization is, forgetting an overall normalization factor: 

g(p T ;Ops,cc,n) = p T e la 'ps , p T < n x a PS ; 

= PT (<L_ r , + £T.\a ' PT>nXO PS ; (1) 

\n ~ Ops' 

where the parameters a, n, and Ops are all positive definite. Fitted to various generator level distri- 
butions, this parametrisation shows a nice behaviour in the range in < pj < 50 GeV. Over a wider range, 
the agreement slightly deteriorates, due to the fact that the high pj power law with a constant power is a 
crude approximation. Both PDFs and matrix element's power depend on the scale Q of the process, which 
is related to pj. The fit quality could be improved by introducing a running power law, a(Q 2 ), at the cost of 
additional free parameters. 

Figure Q] shows the parametrisation fitted to distributions obtained using the Monte Carlo event generators 
PYTHIaH], MC0NL0 ifxxh . and two versions of RESBOS fllEGjI]: the default computation, and a computation 



including small-* broadening effects. All distributions are obtained at = 14 TeV, Q = Mz and yz = 0. 
The fit quality is good, with % 2 ~ 1 in all cases. 



It is interesting to study the dependence of the Ops parameter as a function of rapidity, as illustrated in 
Figure [2] The standard RESBOS prediction shows falling values of this parameter at higher rapidity, an effect 
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2 A parametrized form for the p T distribution 




Figure 1: Generator level pj distribution, at = 14 TeV, Q = Mz and yz = , as predicted by PYTHIA [(a)] 
MC®NL0 [(b")l the standard RESBDS[(c)] and by the version including low-x broadening effects [(d)] 
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generally expected from the decreasing phase space on one side of the parton shower. The modified version, 
however, shows an increase of this parameter, resulting from the low-* effects. It would be interesting to 
measure this dependence in the current and forthcoming hadron collider data. 




(a) (b) 



Figure 2: Rapidity dependence of Gps, as predicted by the standard RESBOS, and by the version including low-x 
broadening effects [(a){ mass dependence of Gps, as predicted by PYTHIA|(b)1 



Finally, one can study the Q dependence of the Ops parameter. As stated above, it is naively proportional to 
s/N~ g , the number of gluons emitted in the initial state. N g is, on the other hand and according to the Altarelli- 



Parisi evolution equations [ 13l4l5h . proportional to the logarithm of the scale variation between the original 
proton and the hard process, N g ~ log(<2 2 ). Plotting Ops as a function of Q, taken to be the boson invariant 
mass event by event, shows a behaviour following G>s ~ \f\og(Q 2 ) as expected in this simple picture. 



3 pr spectrum from the angular distributions 
3.1 Methodology 

The proposed measurement procedure is suggested by observing that at given transverse momentum and at 
fixed mass, the Z boson angular distribution can be written as the product of the lepton angular distribution in 
the Z rest frame, and a factor relating the lepton angles in the rest frame and laboraty frame. For simplicity, 
we consider the azimuthal angular distribution only: 

**J„ \d*')„ \d^) n 

Above, (j>* is the azimuthal angle in the rest frame, and A<p = 0i — 02 the angular separation in the laboratory 
frame. The transverse momentum distribution can be inferred by noting that the overall Atp distribution is 
the integral of the above over the pj distribution: 

do f ( do\ do , 

dp T (3) 



<iA0 J \A(j) J pr dpj 
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3 p T spectrum from the angular distributions 



In Eq. [2 the first factor on the right hand side has a well defined form, and can be computed perturbatively. 
In the Collins-Soper fram^ lfl6fl . one finds: 

— ) ~ 1 + cos 2 0* + (--- cos 2 e*)A + cos 0* sine* cos 0* Ai + -sin 2 0*cos20* A 2 (4) 
dQ.J pT 2 2 2 

where the coefficients A, can be calculated perturbatively and are functions of the kinematic variables s, y, 
p T . The second factor is purely kinematic. In the simplest case where the system is purely transverse (all 
rapidities are and momenta are purely transverse), the relation between 0* and A(f) takes the following 
form: 



where j8 = ^r, and from which the the derivative in Equation [2] can be computed. According to the 
above, the A<p distribution is directly sensitive to pj\ small values of Aip indicate large pj values. In the 
general case, the polar decay angles complicate the picture significantly, as at finite but modest pj, small 
A0 values are also obtained from forward decays (|tj | S> 0), leading to a small projection of the lepton pair 
opening angle in the transverse plane. In addition, the expressions have to be integrated over the Z boson 
lineshape, so that in practice we have to extract the factors from large Monte Carlo generated event samples. 
For the sake of simplicity, we stick to Eq.[3]and compute the the A(j) distribution at given pj integrating over 
all other relevant variables: 

7? J = / ( 7? ) dmdydri (6) 

For the analysis, we use a sample of 10 6 events, generated with MC0NL0. While this approximation is not 
optimal as the lepton rapidities are also measured, providing additional information which is not exploited 
here, it is sufficient for the purpose of demonstration. Statistical sensitivities discussed here should thus be 
understood as conservative. 

The integrated azimuthal angular A<p distribution prediction for MC@NLO is shown in Figure [3l requiring 
that both decay leptons have a transverse momentum above 20 GeV and a pseudo-rapidity rj smaller than 
2.5, as generic acceptance cuts applied by the LHC experiments ATLAS and CMS. The corresponding A<p 
distribution vs. the transverse momentum distribution of the Z is shown in Figure [4] It can be interpreted as 
a matrix M, relating a given pj distribution to a Atp distribution, via 

A<j>i=Mij-p T j (7) 

where i and j are the numbers of bin of the A(j) and pj distribution. It is required that each row of M is 
normalized to unity. The basic idea is to use the measured A<p distribution to estimate the pj distribution, 
exploiting their relation through the matrix M. The inverse matrix M~ 1 directly relates pj distribution to 
the measured Aip by a simple matrix multiplication. The matrix M has significant off-diagonal entries and a 
non-uniform distribution of values on the diagonal. Hence, the inversion of M and the statistical fluctuations 
in Atp induce large fluctuations in pj and therefore sub-optimal results. 



1 /2j8 2 -cosA0 



''Defined as the gauge boson rest frame which maximizes the projections of the beam momenta on the z-axis. 



3.2 Expected Precision 
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Figure 3: Opening-angle distribution for two Z boson de- Figure 4: Distribution of the opening-angle between the 
cay leptons. two leptons and the pj distribution of the Z boson. 



Therefore it was chosen to adjust the MC predicted pj distribution iteratively to minimize the difference 
in the corresponding A<p MC distribution, which is calculated via Equation Ul and the measured A<p Data . The 
difference is expressed as a # 2 -value, i.e. 

; a MCj 
V ^Data,i- M ij-PT,j 

L a 2 

t u MC,i 

It was assumed that the statistical uncertainty o is purely due to the measurement, as the Monte Carlo based 
observables can be theoretically defined with infinite statistics. The parameters to be varied in Equation [8] 
are all entries of the pj distribution, i.e. in realistic scenarios more than 30. This relatively large number of 
free fitting parameters dramatically hinders the minimization of Equation [U especially when the statistics of 
the measured A(f) distribution is limited. Hence it was chosen to use f(pr,Ops, defined in EquationQ] 
to model the pj distribution, i.e. using only three free parameters G>s, a and n during the % 2 -minimization 
procedure. It has been shown in Section [2] that the parameterization of the pj spectrum via Equation Q] 
provides an adequate description up to statistics of at least 5.10 4 events, which we assume in the present 
analysis. 



3.2 Expected Precision 

As already mentioned in Section [T] a prominent systematic uncertainty in the standard measurement of the 
pj is induced by the uncertainty on the decay lepton momentum measurement. In contrast, the presented 
approach has only a very weak dependence on the momentum measurement, but relies on the measurement 
of A<I>, which has in general extremely high precision in most collider experiments. The expected precision 
of modern detectors, like ATLAS or CMS, is a magnitude smaller than the required binning in the AO distri- 
bution. For this reason, the A<I> distribution is not required to be unfolded, which is an additional advantage 
compared to the standard measurement. 
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3 p T spectrum from the angular distributions 



The uncertainties on the lepton transverse momentum measurement can be parameterized to first order as 

Pt —s- a- pj + Gauss(/n, a) 

where a is a scaie parameter, m is an offset parameter of and a is the width of an additional Gaussian (resolu- 
tion) uncertainty. The scale parameter a is assumed to have an uncertainty of 1%, the resolution parameter a 
is assumed to have an uncertainty of 300 MeV. Moreover, we assume a systematic shift of 2 x 10~ 3 rac? on the 
AO measurement. The latter assumption is rather conservative, keeping in mind the precise AO resolution. 

With this assumptions of systematic uncertainties, we can compare both measurement techniques. Table 
Q] shows the comparison of function parameters a,G,n of Equation Q] and its maximum max for different 
measurement methods and systematic uncertainties. Standard denotes that the parameter values have been 
obtained by directly fitting the predicted transverse momentum distribution. AO denotes that the correspond- 
ing values have been obtained with fitting the opening angle distribution. Ideal labels that a perfect detector 
has been assumed, i.e. with perfectly known resolution, while distorted assumes the stated uncertainties on 
the detector resolution. The given values are based on 500.000 selected Z boson events in a specified leptonic 
decay channel, generated with MC@NLO. 

The parameters in the ideal, standard column are the reference values for the comparisons. The values of 
the AO approach agree within their statistical uncertainties to the reference values. Moreover, the assumed 
systematic uncertainty on the AO measurement has no significant effect of the fit results. This is not the case 
for the distorted standard measurement, where a significant difference compared to the ideal measurement 
can be observed. The statistical precision of the AO-approach is reduced less, by a factor of 2. Such a de- 
crease is expected due to the integration over rapidity and the mass of the Z boson and hence not all statistical 
information used. This can be partially recovered when taken the rapidity information into account during 
the fitting procedure, or when restricting the analysis to a smaller rapidity range, e.g. \y\ < 1 .0. Nevertheless, 
the systematic difference between the standard measurement is worse compared to the statistical uncertainty 
of the AO-approach. 



Parameter 


standard 


standard measurement 


A<t> approach 


AO approach 




(ideal) 


(distorted) 


(ideal) 


(distorted) 


max 


3.89 (0.02) 


4.04 (0.02) 


3.90 (0.04) 


3.90 (0.04) 


a 


3.25 (0.03) 


3.10(0.05) 


3.22 (0.8) 


3.23 (0.09) 


a 


0.90 (0.01) 


0.85 (0.01) 


0.89 (0.02) 


0.91 (0.03) 


n 


3.57 (0.02) 


3.94 (0.04) 


3.58 (0.06) 


3.59 (0.06) 



Table 1: Comparison of function parameters a,G,n and the maximum of the function for different measurement 
methods and systematic uncertainties. The statistical uncertainties are given in parenthesis. 



3.3 Comparison of different Monte Carlo Generators 

As previously discussed, the values of the matrix M are predicted by a Monte Carlo generator. In this sec- 
tion, it will be discussed to which extent the final pj measurement is independent of a specific Monte Carlo 
generator program. In order to test this independence, an attempt was made to predict the pj spectrum of 
one generator using the predicted matrix M of a second generator. Again, it was chosen that the matrix is 



3.4 X-Broadening Effects 
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based on MC@NL0. The p T and A<1> spectra used for the comparison and fitting, respectively, are taken from 
HERWIG, RESBOS, and RESBOS including the x-broadening effect. Each spectra is based on 50.000 recon- 
stucted Z boson events, corresponding roughly to an integrated luminosity of S£ = J 100pb _1 at the LHC 
at a collision energy of 7 TeV. The resulting pj spectra, including the functional description obtained with a 
direct fit and via the AO measurement are shown for the different generators in Figure [5] The corresponding 
fitting parameters are shown in Table [2] 



Parameter 


MC0NL0 


HERWIG 


RESBOS 


RESBOS 
(small x) 


GTruth 
Oa* 


3.31 (0.09) 
3.14(0.30) 


3.80 (0.17) 
3.70 (0.44) 


4.27 (0.11) 
4.45 (0.29) 


5.50 (0.13) 
5.70 (0.28) 


®Truth 
(Xa® 


0.90 (0.03) 
0.83 (0.09) 


0.92 (0.04) 
0.88(0.11) 


0.97 (0.03) 
0.99 (0.07) 


1.07 (0.03) 
1.07 (0.07) 


nTruth 
«A* 


3.55 (0.09) 
3.86 (0.20) 


3.81 (0.14) 
4.12(0.39) 


3.62 (0.14) 
3.96 (0.40) 


4.04 (0.25) 
4.84 (0.77) 


maXTruth 
OTflXA<j3 


3.95 (0.07) 
4.11 (0.14) 


4.34 (0.10) 
4.50 (0.22) 


4.48 (0.09) 
4.70 (0.16) 


5.50 (0.13) 
5.78 (0.20) 



Table 2: Comparison of function parameters a,G,n and the maximum of the function for different Monte Carlo 
generator programs. The subscript truth denotes that the function has been directly fitted to the predicted transverse 
momentum distribution, while the subscript A<t> denotes that the corresponding values have been obtained with fitting 
the opening angle distribution. The values correspond to 50.000 selected events. 

The A<I> fitted values agree within their statistical uncertainty to the values, obtained by a direct fit to the 
truth pj distribution. As already mentioned in the previous section, the statistical precision and also the 
systematic differences are mainly due to the integration over the Z boson rapidity. 



3.4 X-Broadening Effects 

As a final example we want to demonstrate that the presented A<I> based pj spectrum measurement can be 
also used to test the x-broadening prediction at LHC for early data, i.e. small integrated luminosities. 

Instead of directly measuring the pj spectrum to test the x-broadening effect, we propose to measure the 
maximum of the pj spectrum for different intervals of the Z-boson rapidity. Figure [6] shows the predicted pj 
spectra for different yz- Larger Z boson rapidities yz test smaller x-regimes of the interacting partons. Hence 
it is expected that the x-broadening enhances for larger yz values, i.e. the maximum the pj shifts to larger 
values. The measurement of the maximum dependence of the pj on yz does not only allow to see a possible 
x-broadening effect, but also to constrain some model parameters. 

To test the A<J> based measurement, we again assume a statistics of 50.000 reconstructed Z boson decays, 
distributed in five rapitity intervals ([0,0.5], [0.5,1.0], [1.0,1.5], [1.5,2.0] and [2.0,2.5]). In each interval, 
we perform the A<t> based fit and extract the maximum of the corresponding pj distribution. The results 
are shown in Figure |7] It becomes evident that we can distinguish between the standard prediction and the 
small-x prediction using the AO based fit on a relatively small data sample. 
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3 p T spectrum from the angular distributions 





Figure 5: Fitted function for different Monte Carlo generator programs: MC0NL0 (a), HERWIG (b),RESB0S (c), RESBOS 
with x-broading (d). The dashed function has been directly fitted to the Monte Carlo predicted transverse momentum 
distribution, while the solid function line has been obtained with by the opening angle distribution. The errors and 
therefore the quoted likelihoods correspond to 50.000 reconstructed events.. 
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Figure 6: Transverse Momentum Spectrum for different Figure 7: X-Broadening determination with 50.000 recon- 
rapidity ranges of the Z boson structed events. 



4 Perspectives 

The present work has two components. First, a phenomenological, three-parameter parametrization of the 
heavy boson pj distribution was introduced, which proved sufficiently versatile to describe the available 
theoretical predictions. This function has many applications, in quantifying differences between predictions 
and assisting the measurement procedure. 

Secondly, we propose a measurement method that is free of any energy resolution systematics, while remain- 
ing a measurement of the pj distribution. The result can thus be directly compared to theoretical predictions 
of this quantity. We presented here a simplified version of this algorithm, based on a map representing the A(j) 
distribution at given pj. The map was integrated over rapidity for simplicity of this presentation, at the cost 
of some statistical power. A complete treatment will have to account for the lepton pseudo-rapidity event by 
event. 

The computation of this map was performed by Monte Carlo simulations, and relies on the well known Z 
boson mass distribution, on kinematics, and on the 0* distribution in the Collins-Soper frame, which can be 
computed perturbatively with good precision. The uncertainty induced by this assumption can be expected 
to be small, but will have to be quantified by further study. A realistic measurement will need sufficient 
statistics, typically 0(50000) events, and would benefit from an analytical calculation of the map. 

In summary, we have proposed here a method that takes as input reliably computed quantities on the the- 
oretical side, and precisely measured angles on the experimental side. The resolution of the lepton energy 
measurment enters only through the kinematic selections, with negligible effect. The statistical power of 
the method is about a factor two less than a direct measurement, but we expect that the reduced systematic 
uncertainties involved here will compensate for this in the long term. 
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